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Abstract. We consider the defect production of a quantum system, initially prepared 
in a current-carrying non-equilibrium state, during its unitary driving through a 
quantum critical point. At low values of the initial current, the quantum Kibble- 
Zurek scaling for the production of defects is recovered. However, at large values of 
the initial current, i.e., very far from an initial equilibrium situation, a universal scaling 
of the defect production is obtained which shows an algebraic dependence with respect 
to the initial current value. These scaling predictions are demonstrated by the exactly 
solvable Ising quantum chain where the current-carrying state is selected through the 
imposition of a Dzyaloshinskii-Moriya interaction term. 


1. Introduction 

When a system is driven adiabatically close to a critical point or through a gapless 
phase, the divergence of its intrinsic relaxation time leads to a complete breakdown 
of the adiabatic condition, no matter how slow the driving is. In the vicinity of the 
gapless point the dynamics switches from an adiabatic to a sudden quench regime. 
Consequently, as we drive the system from its initial ground state closer and closer to 
the gapless point, the departure from the adiabatic evolution induces transitions towards 
excited states which are seen as a proliferation of topological defects. This proliferation 
ultimately generates a hnal state which differs signihcantly from the naively expected 
ground state associated to the hnal Hamiltonian. For a slow driving rate, generalizing 
to the quantum situation the classical Kibble-Zurek Mechanism (KZM) [1] , the density 
of such defects is expected to be given by a power-law function of the driving rate with 
an exponent related to the quantum critical point exponents laEllllElE]. 

The scaling argument goes as follows. For a linear ramping through an isolated 
quantum critical point controlled by the parameter e{t) ~ t/rq, the system, after 
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following adiabatically the ramping far away from the critical point, will suddenly freeze 
out when it gets close enough to the critical locus. This happens at a typical time scale 
tkz which is deduced self-consistently by equating the intrinsic relaxation time with 
the inverse of the instantaneous energy gap A(t) at t = tkz, where the typical scaling 
of the inverse gap, /S.~^{tkz) ~ \^ijKz)\~^'') is set by the deviation from the critical 
point e{TKz) at that time. One obtains tkz ^ accordingly the associated 

typical length scale fxz ~ ^kz ~ which gives an estimate of the defect density 

as Hexc ~ iffz ~ ^^dv/{i+zu)^ This scaling prediction, based on an adiabatic-sudden- 
adiabatic evolution scenario [3], has been tested numerically and by analytical means 
in a great variety of models; see [7[ IH] for extensive reviews. In particular, in a number 
of integrable models, the KZM can be derived exactly from the mapping to a set 
of independent two-level systems, each of them undergoing Landau-Zener (LZ) anti¬ 
crossings [9]. Integrating the LZ transition probabilities over all modes leads finally 
to the KZM prediction for the density of excitations 0171 HD]. Initially introduced for 
homogeneous systems, the KZM has since been generalized to inhomogeneous situations 
such as those generated by the release of a power-law conhning potential m or by the 
propagation of a domain wall or critical front ini. The KZM has also been used in 
spatially inhomogeneous situations to describe symmetry-breaking phase transitions in 
space [13]. 

In all these cases the non-equilibrium situation generated by the temporal variation 
of a Hamiltonian parameter is obtained from an initial equilibrium state (generally the 
ground state since the system is supposed to be at zero temperature). However, many 
situations of physical and technological interest have to do with starting states that 
are intrinsically out-of-equilibrium, for example, situations where there is a macroscopic 
current flowing through the system as a result of a coupling with two different baths or 
reservoirs. Starting the driving process from such an excited state will strongly affect 
the way the defects are generated especially when the dynamical system comes close to 
a gapless point. If the initial current density is small enough, one expects to recover a 
density of defects which is governed by (almost) the usual KZM prediction. However, 
it is possible that this prediction could break down completely at large initial current 
values. It is the aim of this paper to clarify the issue. 

In a one-dimensional system, a current-carrying state can be prepared by the 
contact of the system at its boundaries with reservoirs (or heat baths) at different 
chemical potentials (or temperatures). In such a case, the system steady state will 
in general be a statistical mixture, well described close to equilibrium (small gradient 
of chemical potential or temperature) by a MacLennan-Zubarev density matrix [T3]. 
Far away from equilibrium, there is no general prediction for the density matrix of 
such a steady state although some exact results in terms of Matrix Product States 
have recently been obtained for integrable models [I5l [I6j. Another generic situation 
where such a current-carrying state emerges asymptotically in time is the case of the 
relaxation of a one-dimensional system in which initially the left and right halves have 
been set to different temperatures, or chemical potentials, and then glued together 
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by a local coupling. It was shown on quasi-free systems (such as the quantum 
XY spin chain) that the steady state reached from this type of initial state by a 
unitary evolution is effectively described by a generalized Gibbs distribution of the form 
Q-h{H-\Y) pj jg ^pg Hamiltonian of the chain, /3 = (/di + /3i^)/2 the average inverse 

temperature, A = — (dR)//d the non-equilibrium driving force, and Y a long-ranged 

conjugated current operator (which describes an inhnite set of conserved quantities) 
[HI- Moreover, in one-dimensional critical systems it has been proven by conformal held 
theory techniques that the steady-state currents are universal, depending only on the 
central charges of the theories and on the external temperature or chemical potential bias 
[IS] . In the zero temperature limit the generalized Gibbs or MacLennan-Zubarev density 
matrix p ~ reduces to the ground state projector |GS'y)(GS'y| associated 

to the effective Hamiltonian H — XY. This is basically the physical justihcation for 
the use of the Lagrange multiplier method developed in [T9l EO]. Specihcally, to the 
Hamiltonian H of the considered system one adds a term — AJ proportional to a given 
current operator J (associated to a given conserved quantity). The ground state of the 
effective Hamiltonian H — AJ will then carry a non-vanishing mean current (J) 7 ^ 0 
at sufficiently large values of A. This effective ground state is then interpreted as a 
non-equilibrium, current-carrying, state of the original model described by H. Such 
an effective approach is believed to capture, at least locally, the essential features of a 
system coupled to two different quantum reservoirs at sufficiently low temperature. In 
fact, it is known to yield an exact description of steady states for energy transport in 
critical systems, where J becomes the total momentum operator [T 8 | 1^ . This latter 
observation suggests that analysis of JT — AJ is particularly relevant for studying scaling 
exponents. 

In the following we will apply this strategy to the Ising quantum chain to select 
a state carrying an energy density current. Given that state, we will drive the Ising 
chain through its quantum critical point and focus on the asymptotic defect generation. 
Finally, after the presentation of the model and its explicit solution, we will give an LZ 
calculation leading to the density of defects and extract from it a general argument for 
the generation of defects in such current-carrying situations. 

2. Model 

We start by considering the transverse Ising Hamiltonian on a one-dimensional lattice 
of L sites with periodic boundary conditions: 

( 1 ) 

I I 

Here the afs are the usual Pauli spin matrices and h is a held favouring alignment in 
the 2 ; direction. For dehniteness we take h > 0. It is straightforward to show that an 
energy current in this model can be dehned as 

i i 


( 2 ) 
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with current conservation reflected by [iJ, J] = 0. Following the approach outlined in 
the introduction, we argue that a current-carrying excited state of H can be generated 
as the ground state of the effective Hamiltonian 

Hj = H -\J (3) 

where the A J term introduces an interaction of Dzyaloshinskii-Moriya form and, without 
loss of generality, we assume A > 0. Such an effective Hamiltonian was previously 
studied in, e.g., [191 l22l |23] and, as we recap below, can be brought into diagonal free- 
fermion form by a series of exact transformationsj|] 

Firstly, using the standard Jordan-Wigner transformation 


af = c\ exp ^ c]c^ 


a I =exp 


Cl 


(4) 

(5) 


the complete Hamiltonian can be written as 

Hj= [( 44+1 + 4ci+i - q 4 +i - 

i 

+ \ (Scjc, - l) + ^ (Iq+i + Qc|+i) 


( 6 ) 


where the q’s are fermion operators. The next step is a Fourier transform to wave 
fermions 

1 


H — 


Cl = 


Vl 


Vl 


E 

k 


ate 


ikl 


] ^-ikl 


a'^e 


which, after some manipulation using fermion anti-commutation rules, yields 


Hj = 




sin k {h -\-2 cos k)] a\ak + i sm.k{a\Q^_^ + akcy.-k) + 


h 


Finally, we perform a Bogoliubov-type similarity transform 
ak =r]k cos Uk - ir]l_^ sin Uk 
alfc = -%sina;fc -h coscafc- 


(7) 

( 8 ) 

(9) 


( 10 ) 

( 11 ) 


Here Uk is assumed odd in k, i.e., u^k = —^k which implies that terms of the form 

rj-kVkhXsink x i coscufc sincjfc (and conjugate) are also odd and cancel in the sum. The 

off-diagonal terms in the remaining part are cancelled by the choice 

2 sin /c , , 

tan2a;fc = ——--. (12) 

n. -|- 2 cos k 


f Subtleties relating to the boundary conditions under these transformations are not expected to be 
relevant in the thermodynamic limit; see, e.g., |24] for more detailed analysis. 
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A 

Figure 1. Phase diagram corresponding to Hamiltonian Hj ([^. The mean current 
is non-zero for A > Xc{h) = max(l, 2/h). In later sections we use this construction to 
consider a quench starting from a current-carrying initial state with h > 2. 


Note that this is exactly the same Bogoliubov angle as in the A = 0 case which is, in 
fact, to be expected since the cnrrent operator commntes with the Hamiltonian. With 
the natnral choice that 2uk is in the same qnadrant as the point {h + 2 cos fc, 2 sin k) so 
that ujk takes the sign of k, the Hamiltonian Hj can then be written in diagonal form 
as 

Hj = ^ rjlrjk f hX sin k -|- \J{h-\- 2 cos kY -\- A sin^ k j -|- const. (13) 

We observe immediately that A 7 ^ 0 breaks the k ^ —k symmetry of the spectrnm; 
see, e.g., [T^]. Following |22], we show in £gnre[^the resnlting phase-diagram in h — A 
space noting the slightly different parameterization of onr Hamiltonian to that in the 
literatnre. As indicated by the spectrnm in figure for A > \c{h) = max(l, 2 /h) the 
ground state of Hj consists of a band of Filed current-carrying modes between k- and 
/c+ where 


, -2± W(A2h2 -4)(A2 - 1) 

cos/c± =-(14) 

In concluding this section, we emphasize that the dynamics is unaltered by the 
addition of the Dzyaloshinskii-Moriya term (since A J commutes with H) but the effective 
Hamiltonian Hj (with its A-dependent ground state) can be used as a tool to generate a 
current-carrying initial state. Our programme in the following section is to analyse the 
creation of defects when the system starts from such a current-carrying initial state and 
is then quenched across the critical line which, with our parameterization, is at h = 2 . A 
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k 

Figure 2. Spectrum = h\ sin k + ^{h + 2 cos + 4 sin^ k for /i = 3, A = 2 (solid 
red line) with zero energy line (dashed green) shown for comparison. States between 
k- and k-^. are filled current-carrying modes. 


related discussion on the properties of the entanglement entropy under a similar quench 
can be found in [23]. In fact, there is another critical line at h = —2 but as we restrict 
ourselves throughout to h > 0, our quench protocol never crosses this. 


3. Calculation of defect production 


3.1. Details of dynamics in Heisenberg picture 


In order to analyse the defect production it is helpful to follow the seminal paper of 
Barouch, McCoy and Dresden 1271 and consider the dynamics from the Heisenberg 
viewpoint. We start by writing 0 as a sum over positive modes p, 

L/2 

a = E (15) 

p=l 


where 


Hp = {h -\- 2 cos k) [alcik + cxlfeCt-fe] + hX sin k[alak — 

2i sin k[al,a^_f^-i-aka^k] + h (16) 

with k = {2tt/L)p. With the obvious choice of basis {|0), a|,al^|0), Q!|,|0), al^lO)} 
in which |0) is the vacuum state of the a fermions, we then have the 4x4 matrix 
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representation 

^ h 2i sin k 0 

^ —2ismk 4cosfc + 3/i 0 

^0 0 hX sin fc + 2 cos k + 2h 

^00 0 


0 

0 

0 


\ 


.(17) 


—/lA sin fc + 2 COS fc + 2/i y 


We note that the current-carrying states a|,|0) and «lfc|0) are completely decoupled from 
the other states and the diagonal structure of their submatrix indicates the conservation 
of current for constant held. 

The time evolution matrix, Up{t), in the Heisenberg picture obeys the {h = 1) 
equation 




(18) 


with initial condition Up{to) = I. Here I is the 4x4 identity matrix and we have 
explicitly now included the time-dependence in Hp{t) to allow for the time-dependent 
h{t) which will be of interest in the following. After the system has reached a current- 
carrying steady state corresponding to a non-zero A, we consider a quench in h{t) with 
the usual A = 0 dynamics so that the non-trivial part of (18) reduces to an equation in 
the 2x2 basis {|0), |2) = al,alfc|0)}: 


d_ / Unit) Unit) \ Unit) Unit) ^ / Ht) 2ismk \ 

dt y U2iit) U22it) j y U2iit) U22it) j Y -2ismk 4cos -F 3h(f) y 

where we have suppressed the p subscript in the matrix elements for notational brevity. 
This system of coupled hrst-order differential equations, easily yields decoupled second- 
order ones. For example, we have 


it/" = h'Un + hU[^ - i2i sin k)U [2 

= h'Un + hU'^i f (2i sin k) [(2i sin k)Un + (4 cos k 3h)Ui2] 

- hUn 


= h'Un + hU[i — (4isin^ k)Un — 2sin/c(4cosfc -|- 3h) 


—2i sin k 


= (4 cos k -|- 4h)t/(^ -|- [h' — 4i sin^ k + ii^ cos k -\- 3h)li\Un 
with initial condition Unito) = 1 and = —ihito). Similarly, one hnds 

iU'f 2 = (4 cos k -|- 4:h)U'^2 + — 4i sin^ fc -|- i(4 cos k -|- 3h)h]Ui2 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


with initial condition Unito) = 0, ^( 2 (^ 0 ) = 2sinfc. The differential equations for U 21 
and U 22 are identical to those for Un and Un respectively but with different initial 
conditions: t/ 2 i(to) = 0, U^iito) = —2smk, U 22 ito) = 1, ^ 22 (^ 0 ) = —f[4cosfc-|-3h(fo)]- 
For certain choices of h(f) the corresponding differential equations can be solved 
analytically (at least with the aid of a suitable computer algebra package) in terms 
of extremely tortuous combinations of hypergeometric/special functions. However, our 
focus here is rather on using properties of the solutions to extract the scaling of the 
defect production when the system is quenched across the critical line. Specihcally, we 
hrst calculate the defect density starting from an initial state with A below the critical 
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value (i.e., a steady state with zero current) and then demonstrate how the apparently 
simple change in the analysis required for a current-carrying initial state can lead to a 
dramatic change in the results. 

We consider quenching the system from above the critical line {h = 2) at some 
initial time t^ to below the critical line at some hnal time tf, according to a given smooth 
protocol (with the shorthand dehnitions ho := h{to) and hf := h{tf) now introduced). 
For A < 1, since there is no current, the system starts in the ground state of the p 
fermions which we denote as 


W'(«o)) = ni°‘(*»))- (25) 

k 

where the |Ofc(fo)) are the vacuum states associated to the diagonal fermions at time to 
such that r]k{ho)\Ok{to)) = 0. 

To calculate the production of defects we need to consider the expectation of 
pl{hf)rik{hf) with respect to the time-evolved ground-state |'0(t/)) = U^{tf)\'ijj{to)). 
To understand this, recall that at any time t the system has a held value h{t) and is 
diagonalized in terms of the operators r]l{h(t))rik{h(t)). The associated (adiabatically 
expected) ground state is the vacuum state with respect to these fermions (since the 
single particle spectrum is positive). As a consequence, the number of defects is just 
the number of fermions on top of the instantaneous vacuum. At very low helds, 
h ~ 0, the total number of defects Ylkhki^)hk{h) reduces to the kink number operator 
|^j(l — see [TJ HI] for more explanation. Since the dynamics is most easily 

expressed in terms of the a fermions we use the time-dependent version of the inverse 
Bogoliubov transformation to write 

vlihf)vk{hf) = -is{hf)c{hf)a-kak+s{hffa-kalk + c{hfYalak+is{hf)c{hf)alal^{ 26 ) 


where s{hf) and c{hf) denote trigonometric functions (sine and cosine respectively) of 
the Bogoliubov angle evaluated at held hf. In the remainder of this subsection we 
suppress all hf and tf arguments and indicate explicitly only the initial-time quantities. 
Now, working in the 2x2 basis introduced above, we have 


UriUtU' = 


Un (7i2 


U 2 I U 22 



TJ* 

'^21 

TJ* TJ* 
'^12 '^22 


(27) 


Using the Bogoliubov transformation at to the expression in (27) can then be re¬ 
written in terms of the operators Pkiho) and //^(ho). It turns out that the only terms 
coupling |0(fo)) and (0(fo)| are those proportional to 7)-k{ho)'rf_j^{ho) and the result is 

{^{to)\Ur,lpuU^\^{to)) = c(ho)2[s2|UnP + iscU^^Ul^ - iscUl^U^^ + c^\Ui2?] 

+ *'S(ho)c(ho) [ 5 ^ 7 / 211/11 -|- iscUiiU22 — iscTJi2TJ2i + (?Ui2U22\ 

- is{ho)c{ho)[s^Ul^U 2 l + ISCUI 2 U 2 I - ISCUI^U 22 + C^Ul 2 U 22 ] 

+ s{hof[s^\U2l\^ + lScU2lU;2 - IScU;^U22 + C^\U22?] (28) 

= \c{ho){sUii - icUi 2 ) - is{ho){sU 2 i - icU 22 )\‘^- (29) 
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Figure 3. Numerical evaluation of (31) for quench protocol h[t) = 2 — t/rq from 
to = —tq to t/ = tq with TQ = 10 (red + symbols) and tq = 100 (green x symbols). 
The horizontal dashed lines show the small-A scaling prediction of (37), the diagonal 
dashed line is the large-A prediction of (40); the crossover between the regimes is 
well-described by the error function (solid lines). 


If Af defects are generated during the quench then, in the thermodynamic limit where 
the sum over modes becomes an integral, the defect density is hnally given by 


nexc = lim ^ - 

L —^oo J_j TT 


\c{hQ){sUii-icUi2)-is{hQ){sU2i-icU22)\^ cifc.(30) 


Crucially, for A > 1 the only difference is the initial state - recall that the dynamics 
is unchanged. If we start in the current-carrying regime, then states between k- and 
are occupied and since the dynamics of these current-carrying modes is decoupled, 
we argue that defect production (corresponding to the production of fermion pairs with 
equal and opposite momenta) can only occur for momenta outside this range. Repeating 
the calculations leading to (30), the only difference is a change in the limits of integration 
so that the result is replaced by 


1 

nexc = - / \c{ho){sUu - icUi 2 ) - is{ho){sU 2 i - icU 22 )\‘^ dk 

^ Jo 

+ - f \c{ho){sUu - icUu) - is{hQ){sU 2 i - icU 22 )\‘^ dk. (31) 
^ J-k. 

This integral can be evaluated numerically (see figure using the solutions of the 
differential equations for the matrix elements of U and, signihcantly, the scaling form 
can also be predicted by an LZ argument as shown in the next subsection. 
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3.2. Mapping to a set of Landau-Zener transitions 


If one restricts the matrices Hp (17) associated to the Ising Hamiltonian H to the non¬ 


trivial sector {|0), |2)} the system maps to a set of independent two-level systems each 
described by the Hamiltonian 

H{k, t) = [2 cos k -\- 2h(t)] I — [e{t) -f b{k)] -\- A{k) (32) 

where the a’s are Panli matrices as before, I is here the 2x2 identity matrix, and the 
coefficients are given by e{t) = h{t) — 2, h{k) = 4cos^(/c/2), and A(/c) = —2sinfc. The 
instantaneous eigenvalues are [2 cos k -|- 2h{t)] ± \/[e{t) -|- b{k)Y + associated to 

the instantaneous eigenvectors | ± {t)). The dephasing factor h{k) can be absorbed by 
a redefinition of a local time, t —)■ tk, for each mode k [7]. 

For a driving e{t) = —t/rq starting deep in the disordered phase (h S> 2) and ending 
deep in the ordered phase (h ~ 0), the main contribution to the excitation density comes 
from the modes close to the Fermi point kp = with an excitation probability mi 

(33) 


— g-7rTQA2(A:) _ ^-i-jTTQsm^k ^ g-47rrQ|fc-7rp 


This excitation probability is substantial only for those modes where TqA'^^k) -C 1, that 
is, in a region around the Fermi point of size \k — kp\ ~ which shrinks towards kp as 
the ramping gets slower. Since the initial state from which we start the ramping carries a 
non-vanishing current generated by the population of the negative modes within a region 
[k_,kp] of the first Brillouin zone, the modes k G [/c_,fc+] are dynamically protected 


thanks to the diagonal dynamics (17). No excitation pairs with momenta —fc}, 

where k G [fc_,fc+], can be generated in the course of time. Consequently, the defect 
density from the current-carrying initial state is given by 

1 


^exc 

TT 


r-k+ i>TT 

/ Pkdk+ - 

'0 ^ J-k. 


Pk dk . 


(34) 


Since our initial state has a very large value of the transverse held h, the boundary mode 
—/c_i_ is always far away from the Fermi point kp = tt and, since the main contribution 

_ 1 /Q 

to Pk comes from the region \k — kp\ ^ Tq ' , one can simply omit the hrst integral and 
write 

^ ^ Pkdk . (35) 


n ~ — 

"^exc — 

TT 


|2 


Plugging Pk ^ e in the previous equation, one hnally obtains the excitation 

density 


T? 


, 1/2 


er: 


t[y4irTQ(ir + !_)] , 


(36) 


where erf( 2 ;) is the error function. Analysis of the error function presents two limiting 
cases: 


If (Stttq) ^a tt -|- , then we have 


J- -1/2 

nr\ _^7- ' 


( 37 ) 
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in agreement with the standard Kibble-Zurek result (remembering here that u = 
z = d = 1). Note that, in this limit of slow quenching (large tq, small A), there is 
no dependence on A. 

If (Stttq)"^/^ vr + /c_, then 


T? ~ 


47rr, 


1/2 


Q 


2 _ 1 

+ kJ) ~ —(vr + kJ). 

TT TT 


(38) 


Furthermore, small tt + corresponds to A large and expanding (14) in this limit 
gives 


TT 




^ 2 
ho 


(39) 


where ho is the value of the field at the start of the quench. (Note that, since we 
start in the ordered state, ho > 2 and the term in the bracket is guaranteed to be 
positive.) Hence, we finally get 


Hf, 


ttA 


— ^ I 1 — ^ 


he 


(40) 


We see that, in this large A limit, there is no dependence on tq but the initial field 
ho does plays a role since it controls the initial current. Of course, for a quench 
starting far away from the critical line we have ho —t cx) and this term drops out. 


4. Summary and outlook 

To impose a hnite current on the system, we have to populate the vacuum state with 
modes (current carriers) within a finite set Ij such that the new state is given by 
rifce/j hi| 0 )- Now, if the set Ij has no significant overlap with the critical domain 

_1 /o 

\k — kF\ ~ Tq ' , the excitation density n^xc will be given by the KZM prediction 
Tiexc ^ (with v = z = d = 1 for our model). On the contrary, when there is 

a significant overlap between the two domains, which is exactly what happens at high 
currents, the excitation density will be lowered. At very high current values, the scaling 
of the defect production is given by n^xc ~ A“^. The reason is that, in this case, the 
Lagrange multiplier has to be very high, A ||-H||, and the dominant contribution to 
the effective Hamiltonian H — XJ comes from the current term, the Hamiltonian H 
itself being a small perturbationj§] Consequently, the effective spectrum has a dominant 
contribution of the form hX sin k and, in the absence of H, the ground state is given by 
occupying all negative modes k G [—tt, 0] so there is no possibility of exciting the system: 
all modes are protected. Now, when we add the Hamiltonian H itself, it will slightly 
shift the single-particle spectrum by a term h (for h S> = 2) resulting in hX sin k + h. 
Modes close to the Fermi point kp = — vr will be unoccupied up to the point k- where 
hXsmk_ -|- h = 0, that is, to leading order up to = —vr -|- 1/A. We are then left 

§ The work of [TS] suggests that such large A values can indeed be physically relevant, even for low 
temperatures. 
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close to kp = —7T with an unoccupied domain of size \k — kpl = 1/A that can be excited 
during the quench, leading hnally to = y J^_if\Pkdk = ^ ~ iiA- lower 

helds, the same argument (linearizing the dispersion relation close to kp = —vr) will lead 
to k_ ~ —71 + A ^1 — and consequently to ~ 

In summary, we have studied the defect generation when driving an Ising quantum 
chain through its critical point from an initial state that carries a net energy current. The 
current value is imposed via a Lagrange multiplier held. We have shown that low values 
of the energy current do not affect the density of defects that are generated during the 
crossing of the critical point. In contrast, at large enough currents the current carriers 
are dynamically protected leading to a signihcant suppression of the defect generation. 

It would be very interesting to test the inhuence of an initial current of particles 
on the generation of topological defects in models which do not reduce to a set of 
free particles and where, consequently, the generation of defects does not reduce to 
a Landau-Zener problem. A potential candidate to probe that inhuence is the non- 
integrable Bose-Hubbard model with a slow driving from the Mott insulator phase to 
the superhuid phase. In the zero-current situation there are already some analytical 
KZM predictions based on the fact that the d-dimensional Mott insulator to superhuid 
transition belongs to the universality class of the {d l)-dimensional XY spin model 
[7] . The possibility of extending these KZM predictions to the case of a current-carrying 
initial state, with ((i-l-l)-dimensional classical analogue, is currently under consideration. 
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